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Abstract 

It is shown that the Green's function on a finite lattice in arbitrary space dimension can be 
obtained from that of an infinite lattice by means of translation operator. Explicit examples are 
given for one- and two-dimensional lattices. 



I. INTRODUCTION 



Lattice Green's functions (LGF) are widely used in solving discrete (e.g., Schrodinger) 



icit expressions for the LGF 
^. The form of the 



equations. In the infinite lattice (or thermodynamic) limit expli 
in different space dimensions d are known for a long time 1 
LGF in the direct lattice space ( X ) as obtained by the Fourier transform method is 
usually given by a d— fold integral over the reciprocal space {Q) that can be further reduced 
to known functions. This allows, in particular, to obtain the eigenenergy equation in an 
explicit form too (see, e.g. [3] for an example covering d = 1,2,3 of the Heisenberg spin 
model). An extension of this approach to finite lattices has been proposed in [sl. It is 
based on an approximation shown to account sufficiently well even for subtle details of the 
exact solution, known from the Bethe ansatz in the case of a finite spin chain. Control over 
finite N (the number of lattice sites along one coordinate) corrections turns to be essential 
for a proper description of the wave function symmetry and behavior of the eigenenergies. 
The approach has revealed new classes of excitation modes for lattices in higher dimensions 
and has allowed to find their connection to those known from Bethe solutions in ID. The 
further development of the approach was based on the proof of several exact expressions 
for the finite lattice sums and the discrete Fourier transform in one dimension j^. We have 
recently become aware of the paper [7] where some similar exact formulas have been derived 
and note that our paper in [6] was submitted before publication of 7|]. 

In particular, with the discrete Fourier transform (DFT) we reproduce the exact eigenen- 
ergy equations known from Bethe ansatz. The latter, however, can not be applied to higher 
dimensional lattices. In the present paper we generalize the earlier exact formulas to higher 
dimensions. It turns out that the GF at finite N can be obtained from that of the infi- 
nite lattice under the action of a simple function of the translation operator. This result is 
tested on the d = 1 case studied earlier and an example for the d = 2 case is also given. 
Remarkably, the expression for the finite LGF in two dimensions is obtained in a form of 
a single integral, as the function does not reduce to known ones. This form is also efficient 
in obtaining exact asymptotic expansions in N. The approach is rather general and can be 
applied to different lattice types without limitation on space dimension. 
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II. GREEN'S FUNCTION ON A FINITE SQUARE LATTICE 



The form of the LGF studied below 



" ^ ' ^ "iV f-; A - COS (Q.) - COS (Q,) ' 



^ iV ^ 

?Ti=0 ra=0 

^ A^. m ^ A,, n 
X,F = 0,1,..., AT- 1. 

is the simplest one and, at the same time, contains the important A— parameters related 
to the symmetry of the wave function, as explained in earlier works 5|, IgJ. For simplicity we 
also confine the energy parameter A to the region outside the scattering states, A > 2. Let's 
consider a simple generalization of the formula expressing ([1]) in terms of modified Bessel 
functions 

dze-^^ ^ ^ V+x {z) IgN+Y (z) exp {-ipA^ - iqAy) , (2) 

g=— oop=— oo 

that can be trivially extended to higher dimensions. One clearly sees that in the thermody- 
namic limit ^ oo ( G^^^ {X, Y) ) only terms with p = and q = "survive" . Eq. ([2]) 
then suggests a possibility to act with the "translation operators " like exp(piV^) upon 
the Gj^^ {X, Y) to obtain the full expression at finite N. 

It should be taken into account that X and Y are defined as above and that due to its 
symmetry property the modified Bessel function J„ (z) should actually be written as I\n\ (z). 
In this form the action of translation operator on the series is more transparent. Then by 
splitting the series in two parts we find from ([2]) that 



G,, (X, Y) = dze-'^ f2 (p + ^ exp (^q (^-z A, + N 



_d_ 

dY 



I\X\ (z) I\Y\ [z) 



+ I dze-^' Yl £ e^P (p (-^^- + ^Jv) ) ^^P ("^ (~^^y + ^ 

Jo q=l ^ ^ / / \ V 



d 

dY 



p=U q 

oo oo 



Z) I\Y\ (z) 



\X\ [Z)11Y\ 



I £ £ ^^p (-'^- + ) ^^p (-^ (-'^^ + ) ^1^1 ^1^1 • 



Taking the geometric series of translation operators, we obtain 

1 exp (iA, - N^) 



Gs, {X,Y) 



+ 



1 - exp (-iA, + N^) 1 - exp (^A, - N^) 



X 



1 



+ 



exp (^A, - N^) 



Gtr{\X\,\Y\). 



(3) 



1 - exp {-iAy + N^) 1 - exp {lAy - A^^) 
Note that the operators in the square brackets do not annihilate each other since they act 
differently on (|X| , . Expressions ([I]) and mean that the LGF at finite N 

is a linear superposition of the limit LGF translated with multiples of the lattice period. 
Generalization of Eq. ([3]) to a lattice in an arbitrary dimension is obvious. With this general 
representation one can straightforwardly obtain the LGF from the integral representation 
of LGF in the thermodynamic limit, e.g. 

''^'^ exp{ip\X\)exp{iq\Y\) dp dq 



^TDL 



sq 



using the property 



d 



A — cosp — cosg 27r27r' 



(4) 



/ ( N— ] exp (CX) = f (NC) exp (CX) 



and that |X — = N — X (where / is a function and C a parameter). One can apply it 
directly to the known expressions of the G^^^ for different lattices. 
In the = 1 case we have 



r<TDL 



i\X\) 



2tt 



exp{iq\X\) dq exp(— t>|X|) 



(5) 



cosh V — cos q 2tt sinh v 

where the energy parameter A has been replaced by coshf . Note that the last equality holds 
also when the modulus sign on the l.h.s. is omitted. Thus, one finds 



Gi, (X) 



1 



+ 



exp (iA 



N— 



) 



1 - exp (-.A + N^) 1 - exp (zA - N^) 



exp 



-f 1X1 



sinhf 



exp {—vX) 



exp {iA + {X-N) v) 
1 - exp (-iA - Nv) ' 1 - exp {iA - Nv) 



+ 



i.e. the result obtained earlier 



sinh V ' 



(6) 



61]. For the sake of comparison to higher dimensional lattice 
we also give the result of applying the operator to the first equality in ([6]) and leading to 
the same result. 



27r 



exp {iqX) 1 dq 

cosh V — cos q 1 — exp (— i A + iqN) 2tt 
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^ /■^'^ exp(-igX) ^ exp {iA + iqN) dq 
Jq cosh V — cos q 1 — exp (i A + iqN) 2tt 
However in the case of square lattice such a direct apphcation of translation operators 

to Glq^{\X\,\Y\) does not reduce to a known function. Generally, G^j^^ is a double 

hypergeometric Appel function j^, that can be reduced to elliptic integrals in more simple 

cases (as the one considered in the present paper). Then we can use the approach defined 

by ([3]) and discussed below. An equivalent way would be to substitu.te in ([2]) the well known 

integral representations of Bessel function involving exponentials 

277 ^ 

In {at) Im (at) = / —In+m i'^cit COS y) cxp {-i {u - m) y) , 

Jo 

to finally reduce Ggq {X, Y) to a single integral. 

In analogy with the d = 1 case, according to ([3]) and (jlj), we obtain (7 = 2/A) 

r (X Y) ^ ^ V 1 V exp{iQ^X + iQyY) 

m=0 n=0 \-%yj 



(-2) = 77- / exp {z COS y) exp {izy) dy, 
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An 



° ^1 — (7 cos 0) 

exp(X/-) , exp(zA,. + (Ar-X)/ 



(7) 



;i - exp (-^A, + Nf -)) (1 - exp (^A, + Nf ")) 

exp (Yf +) ^ exp (^A, + (iV - F) / +) 



;i - exp (-2A, + Nf +)) (1 - exp {lAy + iV/ +)) 



where 



7 COS( 



/(0)^ln _ 

1 + -v/l — (7cos( 

The structure of this exact expression is surprisingly simple and indicates the possibility of 
extension to other lattices. In particular, ([7]) contains the periodic symmetry of the LGF 
represented by the transformation X —>■ {N — X) and Y —>■ {N — Y) . For X,Y << N and 
7 7^ 1 we find in the thermodynamic limit, — 00 : 

X+Y 

G™' (X, Y) = f f I I ,4>. (8) 

'l-{-rcos4>f \l+v'l-(7cos, 



Eq. ([8]), the "byproduct" of ([7]), gives at the same time a convenient integral representation 
for the hypergeometric function 4F3 - the explicit form of the (X, Y) found in , that 

arises in various other contexts, see, e.g. |9|. In particular, for X = Y = we find the 
known result 



Another important observation is that in the N 00 limit the symmetry related information 
contained in the A -parameters seems to be lost. However, recalling our earlier results, the 
slow A^— convergence of the the "correction" terms in the vicinity of the critical condition 
7 = 1 proves that the symmetry of the corresponding solutions of the Schrodinger equation 
on a two-dimensional lattice is retained even in the thermodynamic limit. 

In conclusion, the main result of the paper is the derivation of new exact expressions for 
the Green's function on a finite two-dimensional lattice. It is shown that these expressions 
are the result of applying a simple function of translation operator to the Green function of 
the infinite lattice. The derivation outlines the way for generalizations to other lattice types 
and different space dimensions. The exact formulas reduce to known functions and solutions 
in particular cases. These expressions allow to obtain asymptotic expansions. The latter 
are straightforward away from the critical parameter region and can be found by standard 
methods in the vicinity of critical values of parameters. The integral formulas can be easily 
adapted to the region of the energy parameter corresponding to the "scattered states", i.e. 
A < 2. These further developments will be considered in other publications. 
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